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ABSTRACT 

We present spectral and timing analysis of NuSTAR observations of the accreting X-ray pulsar 
2RXP J130159.6-635806. The source was serendipitously observed during a campaign focused 
on the gamma-ray binary PSRB1259-63 and was later targeted for a dedicated observation. The 
spectrum has a typical shape for accreting X-ray pulsars, consisting of a simple power law with 
an exponential cutoff starting at ~ 7 keV with a folding energy of Ef a id — 18 keV. There is also 
an indication of the presence of a 6.4 keV iron line in the spectrum at the ~ 3 a significance level. 
NuSTAR measurements of the pulsation period reveal that the pulsar has undergone a strong 
and steady spin-up for the last 20 years. The pulsed fraction is estimated to be ~ 80%, and is 
constant with energy up to 40 keV. The power density spectrum shows a break towards higher 
frequencies relative to the current spin period. This, together with steady persistent luminosity, 
points to a long-term mass accretion rate high enough to bring the pulsar out of spin equilibrium. 

Subject headings: pulsars: individual (2RXP J130159.6-635806) - stars: pulsars - X-rays: binaries 
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1. Introduction 


2RXP J130159.6-635806, first discovered by 
the ROSAT observatory, was later rediscovered 
in hard X-rays by the INTEGRAL /IBIS tele¬ 
scope and designated w i th the name IGR J13020 - 
6359 ( Bird et ahl 20061 : iRevnivtsev et al.l 120061) . 
The first comprehensive analysis of the tempo¬ 
ral and spectral X-ray propert i es o f this source 
was done bv IChernvakova et al.l ( 20051) using data 
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from the ASCA , BeppoSAX, INTEGRAL and 
XMM-Newton observatories. In particular, XMM- 
Newton data showed coherent pulsations with a 
period of around 700 s. Joint spectral analysis 
of XMM-Newton and INTEGRAL data demon¬ 
strated that the spectral shape is very typical for 
accretion-powered X-ray pulsars (namely, an ab¬ 
sorbed power law with a high-energy cut-off). 

Based on 2MASS archival data. IChernvakova et al 


( 20051) proposed that the binary companion to 
2RXP J130159.6-635806 is a Be star at a distance 
of 4-7 kpc. This sugge stion was later confirmed 
by IColeiro et al. (( 2013 ). who reported the pres¬ 
ence of emission lines of He I A2.0594 yim and 
Br(7-4) A2.1663 ^m, which are typical for a Be 
star. The spectral type of the optical counterpart 
was determined to be B0.5Ve. The orbital period 
of the binary remains unknown. 

X-ray pulsars in binary systems with Be com¬ 
panions (BeXRPs) typically manifest themselves 
as transient sources through either Type I (pe¬ 
riodic flares related to the periastron passage) 
or Type II outbursts (powerful ra re transien t 
events), or a combination of both (e.g., lReiafl201ll) . 

2RXP J130159.6-635806 shows several differences 
from a standard transient BeXRP. Specifically, it 
has a relatively low persistent flux, a long pulse 
period, and it does not demonstrate either Type I 
or Type II outbursts. IChcrnvakova et al.l (120051 ). 
however, did report some variability of its X-ray 
flux. 

Therefore, there are substantial reasons to con¬ 
sider 2RXP J130159.6-635806 as a member of the 
subclass of persistent BeXRPs ( Reig fe Rochd 
19991) . So far, only a few members of this relatively 
small category of objects have been studied in de¬ 
tail: 4U 0352+309/X Persei, RX J0146. 9+6121 / LSI+ 61 


scope Array (NuSTAR) ( Harrison et al.l 120131) 
during observations of the gamma-ray binary 
PSRB1259-63, with three data sets taken in 2014 
May-June (Chernyakova et al., in prep.). In one 
observation, 2RXP J130159.6-635806 appears in 
the corner of ~ 13' x 13' NuSTAR’ s field of view 
(FOV), and, in two more, the source is at the ex¬ 
treme edge of the FOV. Despite the large off-axis 
angles, the NuSTAR data were successfully used 
to extract coherent pulsations. This motivated 
the NuSTAR team to trigger an on-axis 30 ks ob¬ 
servation of 2RXP J130159.6-635806 in order to 
obtain high-quality data for spectral and timing 
analysis. Table Q] lists the NuSTAR observations 
used in this work. 

NuSTAR carries two co-aligned identical X-ray 
telescopes operating in wide energy band from 3 
to 79 keV with angular resolution of 18” (FWHM) 
and half-power diameter (HPD) of 58”. Spectral 
resolution of 400 eV (FWHM) at 10 keV is pro¬ 
vided by independent focal planes for each tele¬ 
scope, usualy referred as focal plane module A and 
B (FPMA and FPMB). 

NuSTAR data can have systematic positional 
offsets as high as 10". Prior to extraction, we 
therefore corrected the world coordinate system 
(WCS) of the event files for the four observations 
to match the PSRB1259-63 and 2RXP J130159.6- 
635806 centroid positions based on cataloged co¬ 
ordinates. 


Since the NuSTAR PSF has wide wings (IHarrison et al 
2013tlAn et ahll2014l) . we investigated how the sur- 


RX J0440.9+4431, and RX J1037.5-564 (lHaberl et al 
1998 : iReig fc Rochelll999h . 


In this paper, we present results of a compre¬ 
hensive analysis of the temporal and spectral prop¬ 
erties of 2RXP J130159.6-635806 in a broad energy 
range, finding some properties that are very un¬ 
usual for BeXRPs. All errors are quoted at the 
90% confidence level unless otherwise stated. 

2. Observations 

2RXP J130159.6-635806 was initially serendip- 
itously observed with Nuclear Spectroscopic Tele¬ 


lace brightness of the source changes with radius 
in order to determine regions where the source 
dominates over the background. We found that 
2RXP J130159.6-635806 is well above background 
23Sthin 120" (~ 92% of PSF enclosed energy; see, 
e.g., An et al. 20141 ). and that the background 
count rate can be robustly measured at radii > 
200" from the source. Taking this into consider¬ 
ation, we defined the corresponding extraction 
regions shown in Fig. |T] for the 4th (on-axis) 
observation. The other three observations, for 
which 2RXP J130159.6-635806 is far off-axis, have 
been treated similarly. Following NuSTAR rec¬ 
ommended standard practice, we chose the back¬ 
ground regions to be on the same detector chip as 
the source. 
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Fig. 1.— Exposure-corrected FPMA (left) and FPMB (right) images of 2RXP J130159.6-635806 for the 4th 
observation in 3 — 78 keV band. The images have been smoothed by a Gaussian kernel with 3 pixel width 
(1 pix = 3"). Each image is color-coded in logarithmic scale. The color bar on the right shows the units 
of the images expressed in 10~ 4 cts pix -1 s _1 . The solid green circle (120" radius) and the yellow shape 
denote regions for the source and background extraction, respectively. The dashed green circle demonstrates 
an angular distance of 200" from the target. The position of the nearby bright source PSRB1259-63 is 
indicated. 


3. Timing analysis 


2RXP J130159.6-635806 is a known source 
of coherent X-ray pulsations at a period of 
~ 700 s with an ave rage sp i n-up rate of v ~ 


2 x 10 


-13 


Hz 


a r-l 


( Chcrnvakova et al. 20051). 


We performed timing analysis of the NuSTAR 
data u sing the xronos (epoch folding tool ef- 
search; iLcahv et al.l 1983h after barycentering the 
data with barycorr. For each NuSTAR observa¬ 
tion, the pulse period and its uncertainty were 
c alculated fo llowing t he procedure de scribed in 
( Boldin. Tsvgankov. fc Lutovinovll2013 ). Namely, 
a large number (10 4 ) of source light curves were 
simulated, the pulse period of each one was deter¬ 
mined with efsearch, and the distribution of the 
corresponding pulse periods was constructed. The 
mean value of this distribution and its standard 
deviation were taken as the pulse period and its 
la uncertainty, correspondingly. Table [T] lists pe¬ 
riod results derived from the FPMA and FPMB 
combined light curves. The inset of Fig. [2] shows 


the evolution of the spin period as a function of 
time. 

As seen from Table Q] and Fig. [2] all four NuS¬ 
TAR datasets are suitable for pulsation detec¬ 
tion. It is also quite evident that periods recorded 
over the time span of 50 days are not consistent 
with each other, clearly showing a decrease in the 
period. We utilized the 1st and 4th NuSTAR 
observations, which have the most accurate pe¬ 
riod measurements and also span the full dura¬ 
tion of the NuSTAR coverage, to measure a period 
derivative of P = —0.0154(5) s/day, equivalent to 
-1.78(6) x 10 -7 s s-\ or v ~ 4.3 x 10~ 13 Hz s~ 4 
This is in agreement with the IChernvakova et al. 

( 2005h spin-up measurement of the second inter¬ 
val of their data, after the ‘break’ at MJD ~ 51900 
(z> ~ 4 x 10 -13 Hz s _1 ). This is quite remarkable 
since there is almost a decade between the period 
measurements. 
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Table 1: NuSTAR observations 


Seq. 

Num. 

Obs. 

ID 

Start Time 
[ UTC ] 

Exp. 

[ks] 

Net Count Rate a 
[ counts s 1 ] 

Period 13 

[s] 

1 

30002017004 

2014-05-04 10:01 

33.3 

(1.49 ±0.02) x 10" 1 

(2.94 ±0.03) x 10" 1 

643.68 ±0.02 

2 

30002017008 

2014-06-02 19:21 

26.4 

(2.00 ±0.12) x 10" 2 

(7.41 ±0.20) x 10" 2 

643.64 ±0.30 

3 

30002017010 

2014-06-14 17:21 

29.1 

(7.05 ±0.19) x 10" 2 

(1.62 ±0.03) x 10 _1 

643.14 ±0.18 

4 

30001032002 

2014-06-24 00:06 

31.6 

1.342 ±0.007 

1.239 ±0.007 

642.90 ±0.01 


a Net count rate in 3 — 78 keV band for FPMA and FPMB extracted from a circular region with a radius of 120 ,/ . 
b Measured pulsation period for 2RXP J130159.6-635806. 

Note. —The target for the first three observations was PSRB1259-63, which placed 2RXP J130159.6-635806 offset by 9.55 7 
from the optical axis. The 4th observation was taken with 2RXP J130159.6-635806 on-axis. 


3.1. Pulse period long-term evolution 


740 

730 

720 

w 710 “ 


^ 690 " 

j! 680 
3 

O, 670 
660 
650 
640 


&<D 


i— 1 —i— 1 —i— 1 —r 


NuSTAR b N 
j i i.i«!> 


56780 56800 56820 56840 
J_i_I_i_1_i_I_i_L_ 


49000 50000 51000 52000 53000 54000 55000 56000 57000 
Time. MJD 

Fig. 2.— Evolution of the pulse period as a func¬ 
tion of time. Blue circles show period measure¬ 
ments published bv IChernvakova et al. ( 2005h us¬ 
ing ASCA , BeppoSAX and XMM-Newton data; 
values obtained in this work are shown by black 
crosses and red points for data from XMM-Newton 
and NuSTAR , respectively. Dashed lines represent 
linear fits to the period evolution with two differ¬ 
ent spin-up rates (see text for the details). 


2RXP J130159.6-635806 regularly fell into the 
FOV of various X-ray telescopes thanks to ex¬ 
tensive observational campaigns dedicated to 
PSRB1259-63 which is located only 9.55' away. 
This allows us to investigate the long-term evo¬ 
lution of the pulse period. We analyzed the 


XMM-Newton (IJansen et all 120011) archival data 


from 2007-201 1 usi ng the procedures described by 


Chcrnvakova et al, (12005 ) and Science Analysis 


Software (SAS) version 14.0.0. The list of selected 
XMM-Newton observations with corresponding 
period measurements are shown in Table [2] 

Fig. [5] presents the long-term evolution of the 
period, showing that the 2RXP J130159.6-635806 
neutron star has undergone very strong and steady 
spin-up during the last 20 years. The first avail¬ 
able v alue of spi n period, me asured in 1994, is 
735 s (IChernvakova et al.ll2005r) . The most recent 
measurements by NuSTAR , from 2014, show the 
period to be around 643 s (Table [I]). This means 
that during the last ~ 20 years the spin period de¬ 
creased by ~ 92 s, corresponding to a mean spin- 
up rate of ~ 1.4 x 10 -7 s s -1 . Fig. [2] shows a 
ch ange i n the average spin-up rate, first reported 
bv IChernvakova et al. ( 2005h . 

We approximated the long-term period evolu¬ 
tion with a linear function with one change in 
slope. A fit shows that the break occurred at MJD 
51300 ± 217 (mid 1999) with a spin-up rate be¬ 
fore and after the break of (4.3 ± 2.7) x 10 -8 s 
s -1 and (1.774 ±0.003) x 10~ 7 s s -1 , respectively. 
As seen from the fit parameters, the spin-up rate 
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becomes significantly higher after the break. As 
noted above, the inset of Fig. [2] shows that the 
NuSTAR data points fit the long-term spin-up rate 
with remarkably high precision. 


Similar behavior was observed for the X-ray 


pulsar GX 1+4, which showe d steady spin up for 
more than a decade (see, e.g ., Bildsten et al . 1997t 
iGonzalez-Galan et al. 2012 + However, GX 1+4 
belongs to the subclass of accreting X-ray pul¬ 
sars known as symbiotic X-ray binaries (SyXBs). 
For BeXRPs, persistent sources typically demon¬ 
strate pulse periods that are relatively stable. 
Examples include the populati on of Be systems 


in the Small Magellanic Cloud ()Klus et al.112014 1 


and the well-k nown low luminosity G alactic sys¬ 


tem X Persei ( Lutovinov et al. 20121+ Transient 


BeXRPs show stro ng spin-up during T ype I and 
Type II outbursts ( Bildsten et al.1 1997 1 with sig- 
nificant sp in- down episodes in between (see, e.g., 
Postnov et aHl2015l+ Therefore, 2RXP J130159.6- 


635806 is a unique source among the BeXRPs be¬ 
cause it demonstrates steady and high long-term 
spin-up with a relatively low and stable luminos¬ 
ity. 


3.2. Pulse profile and pulsed fraction 

Pulsar pulse profiles and their evolution with 
luminosity and energy band depend on the geo¬ 
metrical and physical properties of the emitting re¬ 
gions in the vicinity of the neutron star. In Fig.[3l 
the NuSTAR pulse profiles of 2RXP J130159.6- 
635806 are shown in three different energy bands: 
3-10, 10-20 and 20-40 keV. The lower panel shows 
“soft” ((10-20)/(3-10) keV) and “hard” ((20- 
40)/(10-20) keV) hardness ratios. 

At all energies, the pulse profile can roughly be 


Table 2: XMM-Newton observations 


Obs. 

ID 

Start Time 
[ UTC ] 

Exp. 

[ks] 

Period 

M 

0504550501 

2007-07-08 12:01 

14.1 

681.20 + 0.15 

0504550601 

2007-07-16 19:59 

55.3 

680.06 + 0.10 

0504550701 

2007-08-17 08:38 

11.4 

680.87 + 0.20 

0653640401 

2011-01-06 17:37 

19.9 

665.65 + 0.20 

0653640501 

2011-02-02 18:59 

26.3 

666.50 + 0.20 

0653640601 

2011-03-04 05:59 

12.9 

664.95 + 0.20 



o 0.5 1 1.5 2 

Pulse phase 


Fig. 3.— NuSTAR 2RXP J130159.6-635806 pulse 
profiles in three energy bands (3-10, 10-20 and 20- 
40 keV), normalized by the mean flux. The lower 
panel shows two hardness ratios, (10-20 keV)/(3- 
10 keV) and (20-40 keV)/(10-20 keV) in blue and 
red, respectively. The profiles are shown twice for 
clarity. 
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divided into one main peak at phases 0.0-0.5 and 
two smaller peaks at phases 0.5-0.75 and 0.75-1.0. 
The main feature that changes with energy is the 
depth of the minimum at phase 0.0. As shown in 
the lower panel of Fig. [3l while the “hard” hard¬ 
ness ratio is almost constant, the “soft” hardness 
ratio shows a maximum at phase 0.0 due to an 
increase in the depth of the minimum at 3-10 
keV. Such behavior is caused by differences in the 
source spectrum with pulse phase (see Section 4). 

Fig. [4] shows the pulsed fraction as a function 
of energy. The pulsed fraction is defined as PF = 
(Anax A m i n )/(I max -f- Anin ) ■ where A max and / m in 
are the maximum and minimum intensities in the 
pulse profile, respectively. Defined in this way, 
the pulsed fraction is very high (about 80%). The 
alternative way to characterize the pulsed fraction 
is the relative Root Mean Square (RMS), which 
can be calculated using the following equation: 


RMS = 


(^Eili( p i- < p >f) 

<P> 


(1) 


where Pi is the background-corrected count rate 
in a given bin of the pulse profile, < P > is the 
count rate averaged over the pulse period, and N 
is the total number of phase bins in the profile 
(N = 30 in our analysis). The RMS deviation 
obtained in this way reflects the variability of the 
source pulse profile in a manner that is not sen¬ 
sitive to outliers like the narrow features seen in 
the profile of 2RXP J130159.6-635806. Therefore, 
this quantity has a value of around 30% that is 
much lower than the classically determined pulsed 
fraction and also independent of the energy band 
(see Fig. [I]). 


It is interesting to note t hat i n contrast to the 
majo rity of X-ray pulsars (jLutovinov &; Tsvgankov 
20091) . 2RXP J130159.6-635806 does not show an 


increase in the pulsed fraction at higher ener¬ 
gies. Such uncharacteristic behavior was previ¬ 
ously observed fo r another persistent BeXRP 
RXJ0440.9+4431 ( Tsvgankov. Krivonos. fe Lutovinov 
20121) . On the other hand, Fig. [3] shows that the 


pulsed fraction increases somewhat with energy if 
A m i n is defined from the the pulse plateau rather 
than from the pulse minimum. In other words, 
the peak-to-plateau difference slightly grows with 
energy. 


3.3. Power Spectrum 


The observed 20-year strong and steady spin-up 
of 2RXP J130159.6-635806 reveals the existence of 
a long-term accelerating torque, which indicates 
that the binary interaction must lead to regular 
accretion over a decade-long time scale, although 
this could certainly be episodic (e.g., at periastron 
passages). The torque can be transferred by mat¬ 
ter accreted from either the disc around the neu¬ 
tron star or a stellar wind from the optical coun¬ 
terpart. Unfortunately, there is no strong observa¬ 
tional evidence allowing us to distinguish between 
these two different accretion channels. In both sce¬ 
narios, this process is defined mainly by the mass 
accr etion rate and the mag netic field strength (see, 

e-g., 


Ghosh fe Lamb .197' 


Due to the unknown distance of 2RXP J130159.6 
635806, the luminosity and mass accretion rate 
are highly uncertain; the magnetic field is also 
unknown since no cyclotron line is found in the 
energy spectrum (Sect. [4]). However, some quali¬ 
tative conclusions about the interaction between 
the accretion disk and the neutron star magneto¬ 
sphere can be made from the noise power spectrum 
of the X-ray pulsar. 

According to the “perturbation propagation” 
model, stochastic variations of viscous stresses in 
the accretion disc cause variati o ns of the mass 
accre tion rate (L vubarskii Il997t Churazo\yet al 


20011) . This, in turn, results in a specific shape of 


the Power Density Spectrum (PDS) of the emerg¬ 
ing light curve. Namely, it will appear as a power 
law with slope -1 to -1.5 (but the exact value is not 
well established for X-ray pulsars) up to the break 
frequency, which is the highest frequen cy that can 
be g enerated in the accretion disk (ILvubarski i 
Il997l) . 

In the case of a highly magnetized neutron 
star, this maximal frequency is limited by the 
Keplerian frequency at the magnetospheric ra¬ 
dius, above whic h one can expect a cutoff in 
the source PDS ( Rcvnivtsev et al.l 12009 ). If the 
source stays in spin equilibrium (corotation), the 
cutoff frequency will coincide with the spin fre¬ 
quency of the pulsar, whereas, in the case of spin- 
up (increased mass accretion rate), the magne¬ 
tosphere will be squeezed and additional noise 
will be generated at higher frequencies. If the 
mass accretion rate is known (i.e., the distance 
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Fig. 4.— The pulsed fraction (red circles) and 
relative RMS (blue squares) of 2RXP J130159.6- 
635806 obtained with NuSTAR as a function of 
energy. 



Fig. 5.— The noise power spectrum of 
2RXP J130159.6-635806 obtained with the NuS¬ 
TAR data in observation 30001032002 (red his¬ 
togram; the 4th observation). The solid line shows 
a broken power law model with the fitted value 
of the break frequency, 0.0066 Hz (vertical dotted 
line). The power law slope above the break is fixed 
to -2. The position of the pulse frequency (0.0015 
Hz) is shown by the vertical dashed line. 


to the source is known), this property of the 
PDS can be used to estimate t he magnetic field 
stren gt h o f t he neutron star dRevnivtsev et al 


2009: Tsvgankov. Krivonos. fe Lutovinov 2012t 


Doroshenko et al.l 1201411 . The appearance of the 


break in the PDS does not necessarily indicate 
that accretion is from a disk. There are wind- 
accreting sources in spin equilibrium also showing 
a brea k in their P DSs a round the pulse frequency 
( Hoshino fe Takeshimal 1993 1. However, the evo¬ 
lution of the PDS shape aw a function of mass 
accretion rate in such systems is not well studied. 

In Fig. El we show the PDS of 2RXP J130159.6- 
635806 obtained with the NuSTAR data in the 
4th observation after subtracting the pulse profile 
folded with the measured period from the light 
curve. The solid line represents the fitting model 
in the form of a broken power law. The mea¬ 
sured break frequency is 0.0066 Hz (shown by dot¬ 
ted line), and it is clear that it is shifted towards 
higher frequencies relative to the spin frequency in 
this observation (0.0015 Hz; shown by dashed ver¬ 
tical line). The power-la w slope above the brea k 


frequency is fixed at -2 ( Revnivtsev et al. 2009?! . 


The best-fit value of the slope below the break 
is —0.33. Given the steady persistent luminosity, 
we conclude that the spin-up observed during last 
~20 years is caused by a long-term mass accretion 
rate that is high enough to squeeze the magneto¬ 
sphere inside the corotational radius. This finding 
confirms the uniquness of 2RXP J130159.6-635806 
among the other X-ray pulsars in binary systems 
with Be companions. 


4. Spectral analysis 

We used nuproducts , a part of the NuSTAR- 
DAS package, to extract source and background 
spectra and to generate NuSTAR response matrix 
(RMF) and effective area (ARF) files for a point 
source. In our analysis we utilized the most recent 
calibration database (CALDB), version 20150316. 
The extracted spectra were then grouped to have 
more than 20 counts per bin using the grppha tool 
from the HEAsoft 6.15.1 package. We fit t he NuS¬ 
TAR spectra using XSPEC version 12.8.1 (Arnaud 
199611 . 
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4.1. Pulse phase-averaged spectroscopy 


According to measurements from the ASCA 
and XMM-Newton observatories, the spectrum 
of 2RXP J130159.6-635806 is characterized by a 
moderate absorption value -/Vh = (2.48 ± 0.07) x 
10 22 cm~ 2 , which i s sta ble over a dozen years 
( Chernvakova et ahl 2005 ). This value was ob¬ 


tained by approximating the source spectra as a 
power law model modified by interstellar absorp¬ 
tion (wabs model in the XSPEC package). Note 
that it is just slightly higher than the value of the 
interstellar hydroge n absorpti on (1. 7 — 1 . 9) x 10 22 
cm -2 determined bv lDickev fc Lockman ( 199 0;) in 
the direction of 2RXP J130159.6-635806. 

In the first three NuSTAR observations 2RXP J130159 
635806 serendipitously appeared highly offset from 
the optical axis, and at the large off-axis angles 
the absolute flux measurements can be affected 
by inaccurate PSF positioning and subsequent in- 
apropriate weighting in the spectrum extraction 
procedure. This leads to large systematics for 
the extracted spectra. Therefore, we restrict de¬ 
tailed analysis of the source spectrum to the data 
obtained in the 4th (on-axis) observation, which 
provides high-quality data. 

In general, the spectrum of 2RXP J130159.6- 
635806 has a shape which is typical for accret¬ 
ing pulsars in binary systems, showing a hard 
spectrum with an exponential cutoff at high ener¬ 
gies. Fig. [6] (a) presents the phase-averaged spec¬ 
trum approximated with the most suitable spec¬ 
tral model determined below. We initially mod¬ 
eled data with a cutoff power-law model (CUT- 
OFFPL in the XSPEC package) 


> 

0) 

2* 


6 ™ 


AE T e E f< 


( 2 ) 


where T is the photon index, Efoiu is the charac¬ 
teristic energy of the cutoff (e.g., the folding en¬ 
ergy), and A is a normalization. This model was 
modified by interstellar absorption in the form of 
the WABS model. As the working energy range 
of the NuSTAR observatory begins at 3 keV, it 
is not very sensitive to measuring low absorption 
columns. Therefore, in the following analysis, the 
interstellar absorption was fixe d to the va l ue Ah 
1 -x 10 22 cm” 2 measured by Ichern vakova ct al 


(12005 1. Note that we simultaneously fitted spec¬ 
tra obtained by both NuSTAR modules. To take 
into account the uncertainty in their relative cali- 
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Energy, keV 

Fig. 6. Top panel (a): Pulse phase-averaged 
spectrum of 2RXP J130159.6-635806 obtained 
with NuSTAR during the 4th observation. Red 
points and blue crosses show module A and B 
data, respectively. The black line represents the 
best fit by the ‘highecut’ model with an iron emis¬ 
sion line at 6.4 keV. Bottom panels (b,c,d) show 
corresponding residuals from the models ‘cutoff’, 
‘highecut’, and ‘highecut’ with the 6.4 keV line, 
respectively (see text and Table [3]). 





























brations, which may be even more of a concern for 
the observations where the source is highly offset 
from the optical axis, we added a cross-calibration 
constant between the modules. The fitting param¬ 
eters for different data sets are shown in Tabic [3] 


The ‘cutoff’ model approximates the source 
spectrum relatively well with x 2 = 992.18 for 921 
degrees-of-freedom (do/; see Table [3]). Neverthe¬ 
less, a wave-like structure is clearly seen in the 
residual panel Fig. [6] (b). To improve the quality 
of the fit, we applied another continuum model 
in the form of a power-law multiplied by a high- 
energy c utoff ( HIGHECUT m odel in the XSPEC 
package; Iwhite et al. 1983} h This model can be 
written as: 


AE~ r x 


1 (E < E cut ) , . 

e -(E-E cM )/E !old ( E > W 


where E cnt is the energy where the cutoff starts. 
As in the previous case, we fixed the interstel¬ 
lar absorption value. Residuals of modeling the 
source spectrum with ‘highecut’ are presented in 
Fig. [6] (c), and best fit parameters are listed in Ta¬ 
ble [31 The ‘highecut’ model significantly improves 
the quality of the fit (x 2 = 929.72 for 920 dof). 
The high-energy cutoff value E cut is found to be 
6.48 ^q' 22 keV, which is significantly l ower than the 
value of ~ 25 keV reported bv IChernvakova et al 


(12005 1 using simultaneous XMM-Newton and IN¬ 
TEGRAL data. The apparent discrepancy is 
probably due to the lower statistical quality of 
the INTEGRAL data and the gap between the 
energy bands covered by the XMM-Newton and 
INTEGRAL observatories. An additional possible 
problem is that the cutoff energy E cut = 6.48 keV 
is very close to the energy of the iron fluorescent 
line at 6.4 keV. The simplistic ‘highecut’ model 
might hide the presence of the emission line in the 
source spectrum. Observed deviations of the mea¬ 
sured spectrum from the model near this energy 
argue in favor of this possibility. 


To investigate this issue, we added a Gaussian 
emission line to the ‘highecut’ model, fixing its 
energy to 6.4 keV and width to 0.1 keV, allow¬ 
ing its normalization to be a free parameter. This 
resulted in an additional improvement of the fit 
t° X 2 = 916.97 for 919 dof for a normalization of 
(2.31±0.70) x 10 -5 ph cm -2 s -1 . The correspond¬ 
ing equivalent width of the line is EW = 44 / 3 ® 
eV (3a error). We determined the significance 


of the line using the XSPEC script simftest with 
4 x 10 4 trials, and found that presumption against 
the null hypothesis, or no line required by the data 
is 3 x 10 -3 , which corresponds to ~ 3 ct line detec¬ 
tion, assuming a normal distribution. The residu¬ 
als of the ‘highecut’ model with the 6.4 keV iron 
line are shown in Fig. [6](d). The model itself, to¬ 
gether with spectral data points, is shown in Fig. [6] 
(a). 

Considering the moderate spectral resolution 
of the NuSTAR observatory and the low signif¬ 
icance of the detected iron line, we investigated 
the possibility of the presence of an iron fluo¬ 
rescence line in the XMM-Newton spectrum of 
2RXP J130159.6-635806. As mentioned above, 
the X-ray pulsar 2RXP J130159.6-635806 was ob¬ 
served with XMM-Newton many times during pro¬ 
grams studying PSRB1259-63. For our purposes, 
we chose the two observations with the longest 
exposures - ObsID 0092820301 (~41.2 ksec) and 
ObsID 0504550601 (~55.3 ksec). As for the NuS¬ 
TAR observations, we modeled the spectra of 
2RXP J130159.6-635806 including the Gaussian 
line at 6.4 keV. For both of the XMM-Newton ob¬ 
servations, we did not find a significant improve¬ 
ment in the fits when the iron line was added and 
obtained a conservative upper limit (3er) for the 
equivalent width of the iron line of 110 eV, which 
is consistent with the NuSTAR results. 


Despite the fact that the first three NuSTAR 
observations were made at large offset angles, and 
the corresponding statistics are significantly lower 
than for the 4th observation, some useful spectral 
information can still be extracted. For these obser¬ 
vations we used the same ‘highecut’ model. Due 
to low statistics and poor fit we fixed power-law 
slope in the 2nd and 3rd data set at T = 1.37 value 
measured in the 4th observation. As seen from 
Table m where the best fit parameters are listed, 
the principal parameters - power-law slope (T), 
cut-off energy (E cut ), and folding energy (E io i d ) 
- of the first three high-offset observations are in 
general agreement with 4th (on-axis) observation 
modeled with ‘highecut’. 


The estimated 2-10 keV flux of 2RXP J130159.6- 
635806 is about 3 x 10 -11 erg s _1 cm -2 dur¬ 
ing our observations, which is in agreement with 
the value of (2 — 3) x 10~ n erg s _1 cm -2 mea¬ 
sured bv IChernvakova et al .1 ( 20051 ). An absence 
of strong transient activity from 2RXP J130159.6- 


9 
















Table 3: Parameters for the 2RXP J130159.6-635806 phase-averaged spectral analysis based on NuSTAR 
observations. 


Seq. 

Num. 

Model a 

Const b 

Photon 

index 

£cut 

[keV] 

-Ufold 

[keV] 

Fluxfj_ 

FPMA 

•10 keV 

FPMB 

xi (dof) 

1 

2 

3 

4 

HI 

HI 

HI 

CU 

HI 

HI+Fe d 

1.24 ±0.02 

0.61 ±0.04 

0.92 ±0.03 

1.00 ±0.01 

1.00 ±0.01 

1.00 ±0.01 

1 4f) +uu/ 

1.37 (frozen) 
1.37 (frozen) 
1.04 ±0.03 

1.32 ±0.03 

1.37 ±0.04 

7.47 +b '55 

r ? n±0.65 
°-‘ z -1.90 

5-68±g;il 

6.48i8;3 

6.94l°;|° 

17.4311;'* 
15.4915+ 
17.4811(35 
12.62 ±0.48 

16.7818:3 

17.96 ±0.87 

z,lu -0.38 

3.2915(47 

911 +0.08 
z - 11 -0.13 

Q 7Q+0.05 

1 ^-0.10 
q 7Q + 0.05 

1 ^-0.10 

Q 7Q+0-05 

1 y -0.10 

2.68!8(£ 

2.0118)95 

1 qc+o.06 
i>yo -0.10 

Q 7Q+0-05 

1 ^-0.11 
q rrQ + 0.05 

1 y -0.11 

Q 70+0-05 
' ^-0.11 

1.06 (458) 

0.92 (124) 

0.81 (264) 

1.07 (921) 

1.01 (920) 

1.00 (919) 


“The xspec spectral model used. ‘CU’: wabs*CUT0ffpl, ‘HI’: wabs*powerlaw*highecut, ‘HI+Fe’: 

WABS(POWERLAW*HIGHECUT±GAu). 

^Constant factor of the FPMB spectrum relative to FPMA. 

“The absorbed flux in units of 10“ n erg s —1 cm -2 . 

r/ The corresponding 6.4 keV iron line parameters are: a = 0.1 keV (fixed), normalization (2.31 ±0.70) X 10 —5 ph cm -2 1 . and 
line equivalent width EW = 44lg® e V (3cr error). 


635806 on long time scales is confirmed with the 
RX TE/ASM and Swift /BAT instruments. Note 
that Chcrnvakova et ahl (2005*) reported on a flar¬ 
ing episode, but the flux of the flare was rela¬ 
tively low raising by a factor of a few (up to 
~ 10 -1 ° erg s -1 cm -2 ). 

Finally, in order to check for the possible 
presence of a cyclotron absorption line in the 
source spectrum, we modified the best-fit model 
(WABS*(POWERLAW*HIGHECUT+GAu)) by in¬ 
cluding an absorption feature in the form of a 
Lorentzian optical depth profile (CYCLABS model 
in XSPEC; Mihara et al. 1990h . The search pro¬ 
cedu re was performed following the prescription 
ffom lTs y gankov fc Lutovinovl (120051 ). Namely, we 
varied the energy of the line over the range be¬ 
tween 5 and 50 keV with 3-keV steps and the line 
width between 2 and 12 keV with 2-keV steps, 
leaving the line depth as a free parameter. We 
did not find strong evidence for a cyclotron line 
in the spectrum of 2RXP J130159.6-635806 (no 
trials gave a significance higher than ~ 2cr). 


4.2. Pulse phase-resolved spectroscopy 

In order to study the evolution of the source 
spectrum over the pulse period, we performed 
pulse phase-resolved spectroscopy using the data 
from the 4th observation. The period was di- 



Pulse phase 


Fig. 7. Parameters of the best fit model 
(WABS* (POWERLAW*HIGHECUT±GAU)) as a 
function of pulse phase in the 4th (on-axis) 
observation. Top: the black histogram shows the 
pulse profile in the entire energy range (duplicated 
in the other two panels). Middle and bottom, 
respectively: cutoff energy (E cut ) and folding 
energy (-E fo i d ). 
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vided into 10 phase bins with zero phase coin¬ 
ciding with the main minimum in the pulse pro¬ 
file. According to the pulse phase-averaged spec¬ 
tral analysis, the fitting model was chosen in the 
form of (WABS*(POWERLAW*HIGHECUT+GAu)). 
However, due to much lower statistics, the pho¬ 
ton index was fixed at the value from the averaged 
spectrum (F = 1.37). We selected this parameter 
due to its virtual constancy over the pulse in our 
preliminary analysis of the same data. The TVh 
value is not very well constrained by the NuSTAR 
data; however we checked that it is consistent with 
being constant in phase, and we fixed it as well. 

The results are shown in Fig. [7] The aver¬ 
age pulse profile of 2RXP J130159.6-635806 across 
the entire NuSTAR energy range is presented in 
the upper panel. The lower panels demonstrate 
the behavior of the two free spectral parameters: 
the cut-off energy (E cut ) and the folding energy 
(Flfoid) ■ The cut-off energy is quite stable over 
the pulse except at the pulse minimum where 
its value increases approximately by a factor of 
two. The folding energy demonstrates an appar¬ 
ent correlation with the pulse intensity, which is 
probably caused by increasing a spectral hardness 
around the pulse maximum. Such behaviour of 
the spectral parameters over the pulse can ex¬ 
plain the corresponding behavior of the hardness 
ratios constructed from the pulse profiles in dif¬ 
ferent energy bands, previously shown in Fig. [3] 
The (20 —40)/(10 —20) keV hardness ratio demon¬ 
strates correlation with the pulse intensity and the 
(10 —20)/(3—10) keV ratio peaks at the pulse min¬ 
imum right in place where -E C ut shifts from ~ 7 to 
~ 16 keV. Finally, it is necessary to note that the 
observed behavior of spectral parameters with the 
pulse phase can be caused by both physical and 
artificial reasons (in particular, due to limitations 
of available data, the adopted spectral model, and 
which model parameters were fixed). More data 
are required to confidently constrain all the spec¬ 
tral parameters and trace their behavior with the 
pulse phase. 

5. Summary 

We summarize the results of spectral and tim¬ 
ing analysis of serendipitous and dedicated NuS¬ 
TAR observations of the accreting X-ray pulsar 
2RXP J130159.6-635806 in 2014 May-June. 


• The source demonstrates strong pulsations 
with a period of ~ 640 s. The ~ 80% pulsed 
fraction is measured to be constant with en¬ 
ergy up to 40 keV. 

• The pulse profile is virtually independent of 
energy and can roughly be divided into one 
main peak at phases 0.0-0.5 and two smaller 
peaks at phases 0.5-0.75 and 0.75-1.0. The 
only feature that is changing with energy is 
the depth of the main minimum at phase 0. 


The measured period shows a significant 
change over the ~ 50 day time span of the 
NuSTAR observations, which is consistent 
with a spin-up rate of v ~ 4.3 x 10 -13 Hz s . 
This rate is in remarka ble agreement wit h 
measu rements made by Chernyakova ct al 


([20051 ) almost a decade ago. 


Toget her with the results of lChernvakova ct al 
( 2005 1. the XMM-Newton data taken in 
2007 and 2011, and the NuSTAR observa¬ 
tions, we show a long-term spin-up trend of 
the source during the last 20 years. During 
the last 15 years, the source has undergone 
a strong and steady spin-up rate at the level 
of P= (1.774 ± 0.003) x 10“ 7 s s" 1 . 


• The power density spectrum of the source 
shows a clear break at 0.0066 Hz, which 
is higher than the period frequency of 
0.0015 Hz. This fact, together with the 
steady persistent luminosity of the source, 
implies that the spin-up observed during the 
last 20 years is likely caused by a long-term 
mass accretion rate high enough to squeeze 
the magnetosphere inside the corotational 
radius, which makes 2RXP J130159.6-635806 
unique among the other X-ray pulsars in bi¬ 
nary systems with Be companions. 


• The phase-averaged spectrum of the source 
has a typical shape for accreting neutron 
stars in binary systems, in particular, for 
X-ray pulsars, and demonstrates an expo¬ 
nential cutoff at high energies. Our best-fit 
model contains an absorbed power-law with 
r ~ 1.4 modified by a high-energy spectral 
drop with a cut-off energy of E cut ~ 7 keV 
and a folding energy of E [ 0 id — 18 keV. The 
spectrum also shows ~ 3cr evidence for an 
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iron 6.4 keV emission line, with an improve- Chernyakova, M., Lutovinov, A., Rodriguez, J., & 

ment in the fit when it is included. Revnivtsev, M. 2005, MNRAS, 364, 455 


• The observed flux corresponds to an unab¬ 
sorbed luminosity in the range ~ (8 — 26) x 
10 34 erg s 1 , assuming a source distance 
betw een 4 and 7 kpc ( Chernyakova et ah] 
These luminosity values imply that 


2005). 


the source is a member of the subcl ass of 
persi stent low luminosity Be systems (jReig 


201 111 . 


• The phase-resolved spectroscopy shows some 
differences in source spectrum with phase. 
The cut-off energy is very stable over the 
pulse except at zero phase where its value 
increases by a factor of two. The appar¬ 
ent correlation of the folding energy with 
pulse intensity is attributed to the change in 
hardness of the source spectrum with orbital 
phase. 


This research has made use of data obtained 
with NuSTAR , a project led by Caltech, funded 
by NASA and managed by NASA/JPL, and 
has utilized the NUSTARDAS software package, 
jointly developed by the ASDC (Italy) and Cal¬ 
tech (USA). This research has also made use of 
data obtained with XMM-Newton, an ESA sci¬ 
ence mission with instruments and contributions 
directly funded by ESA Member States. AL and 
ST acknowledge support from Russian Science 
Foundation (grant 14-12-01287). 

Facilities: NuSTAR, XMM-Newton. 
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